Using Final State Pseudorapidities to Improve s-channel Resonance Observables at 

the LHC 
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We study the use of final state particle pseudorapidity for measurements of s-channel resonances 
at the LHC. Distinguishing the spin of an s-channel resonance can, in principle, be accomplished 
using angular distributions in the centre-of-mass frame, possibly using a centre-edge asymmetry 
measurement, Ace- In addition, forward-backward asymmetry measurements, Afb, can be used 
to distinguish between models of extra neutral gauge bosons. In this note we show how these 
measurements can be improved by using simple methods based on the pseudorapidity of the final 
state particles and present the expected results for Afb and Ace for several representative models. 
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I. INTRODUCTION 

The startup of the CERN Large Hadron Collider 
(LHC) will allow the exploration of the TeV energy 
regime and the testing of the multitude of proposed the- 
ories of physics beyond the Standard Model. Many of 
these theories predict the existence of massive, neutral s- 
channel resonances P, [3, i, i, i, S 0, i, Bill • For some 
models of new neutral gauge bosons {Z'), the LHC is ex- 
pected to have a discovery reach u pwa rds of 5 TeV with 
100 fb~^ of integrated luminosity [11|. This is a signif- 
icant improvement over the current experimental limits 
for most models, which constrain Z' masses to values 
greater than ~ 1 TeV [H [H [ll. 

If a TeV scale s-channel resonance were discovered, the 
immediate task would be to identify its origins. Many 
observables have been proposed to this end, primarily fo- 
cused on the dilepton channel (e and /i) , which would pro- 
duce the cleanest and most easily measured signal for a 
non-leptophobic [lli,i,|l[| with the ATLAS 
and CMS [iTj detectors. The proposed measurements for 
the dilepton channel are the Z' width, total cross sec- 
tion, forward-backward asymmetry {Aps) p^ . central- 
to-edge rapidity ratio fl9| , and a comprehensive analysis 
of all rapidity regions 'id. 21]. 

There are challenges associated with some of these 
measurements that we argue can be alleviated by us- 
ing the pseudorapidity of the final state fermions. In 
particular, we focus on two measurements: determining 
(or at least constraining) the spin of an s-channel reso- 
nance, and determining the forward-backward asymme- 
try of a Z' . Distinguishing whether the new resonance is 
a scalar, such as an R-parity violating sneutrino [ol, [lOj, 
a spin-2 boson^uch as a KK graviton [1, 0, Q , or a spin- 
1 Z' P, [3, H, H H will be challenging and is typically 
determined through the study of the angular distribu- 
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tion in the centre-of-mass frame of the initial state quark 
and anti-quark (cm.) [1, [H HI [13, 111 . The forward- 
backward asymmetry measurement at the LHC has to 
deal with the ambiguity in defining the forward direc- 
tion due to the inability to unambiguously determine the 
direction of the initial state quark in a symmetric proton- 
proton collision. 

Presently, some solutions exist to deal with these chal- 
lenges. To distinguish the spin of the resonance, a centre- 
edge asymmetry, Ace^ |23| can be defined that is sensi- 
tive to the the angular distribution of the events. The 
centre-edge asymmetry is a simple means of binning the 
events in the central and edge regions of cosO* , the cm. 
scattering angle, which will be weighted differently de- 
pending on the angular distribution. This has the benefit 
of eliminating some of the systematic uncertainties of a 
fit to the angular distribution. However, the Ace observ- 
able still relies on boosting the particle four-momentum 
from the lab frame to the cm. frame. 

The forward-backward ambiguity in a symmetric pp 
collision can be resolved by exploiting the fact that 
the valence quarks have, on average, larger momen- 
tum than the sea anti-quarks. The quark direction can 
then be identified with the boost direction of the dilep- 
ton system. [IB]- Restricting the measurement to those 
events that have a large boost (i.e.: \Yzi\ > 0.8) re- 
duces the misidentification of the initial state quarks 
and antiquarks, resulting in greater than 70% of dilepton 
events being correctly identified as being boosted by the 
quark [26] . Both of these methods have been explored in 
great detail and remain the standard approach used in 
the literature [13, [lH, [13, [l^ . 

Both the Ace and the Aps measurements require 
analysis of the centre-of-mass (cm.) angular distribu- 
tion of the dilepton events - directly for Ace, and when 
tagging forward or backward events in Afb- In this 
note we propose a simpler method of measuring these 
asymmetries without reconstructing the angular distribu- 
tions. Specifically, we exploit the direct measurements of 
the lepton pseudorapidities to calculate the observables, 
rather than using derived quantities that may propagate 
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uncertainties into the result. The proposed methods also 
take advantage of the fact that differences in pseudora- 
pidities are Lorentz invariant quantities, so that all cal- 
culations can be performed using quantities measured in 
the lab frame. 

The point of these methods is not to provide new 
phenomenological insight into the models, but rather to 
demonstrate how the use of final state pseudorapidities 
provides a simpler and cleaner means of obtaining the 
Ace and Aps values. The dimuon signal is very clean 
and error propagation should not be a big issue. The 
real power of this approach will be seen when applied to 
heavy quark final states [l^l- In the following sections, 
we give some calculational details which are followed by a 
description of our approaches to the Centre-Edge Asym- 
metry and the Forward-Backward Asymmetry. We con- 
clude with some final comments. 



II. CALCULATIONAL DETAILS 
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FIG. 1: Normalized parton-level angular distribution of 
spin-0 (black), spin-1 (dark grey/orange) and spin-2 (light 
grey /yellow) bosons decaying to fermions. 



The basic ingredients in our calculations are the cross 
sections, cr{pp — > i? where R = Z' , i> or 

G. The cross section for R — Z' \s described by the 
DrelFYan process with the addition of a [Ulillll. 
Analagous expressions for the spin-0 v and spin-2 gravi- 
ton are given in Ref. and Q, respectively. We com- 

puted the cross sections using Monte-Carlo phase space 
integration with weighted events and imposed kinematic 
cuts to take into account detector acceptances, as de- 
scribed in the following sections. 

In our numerical results we take a = 1/128.9, sin^ 9^ = 
0.231, Mz = 91.188 GeV, and Tz = 2.495 GeV [pi. We 
used the CTEQ6M parton distribution functions [30| and 
included a K-factor to account for NLO QCD corrections 
[3l| . We neglected NNLO corrections, which are not 
numerically important to our results (33 . [33| , and final 
state QED radiation effects, which are potentially impor- 
tant [3^ but require a detailed detector level simulation 
that is beyond the scope of the present analysis. The Z' 
widths only include decays to standard model fermions 
and include NLO QCD and electroweak radiative correc- 
tions 1351 . For the v width, we take Fj> = 1 GeV following 
Ref. M|. Expressions for the G width can be found in 
Ref. EI[2l,[3E[33. 



III. SPIN DISCRIMINATION USING 
CENTRE-EDGE ASYMMETRY, Ace 

The parton level angular distributions, da /d cos 6*, of 
the spin-0, -1, -2 bosons, shown in Fig. [H are distinct 
enough that, in principle, such a measurement would 
uniquely identify the spin [§!, ^2^. However, these dis- 
tributions are not directly accessible due to the convo- 
lution with the parton distributions of the protons, the 
boosting of measured lab frame quantities to the centre- 
of-mass frame, detector limitations and finite statistics. 



all of which will make the measurement challenging [22] ■ 
The centre-edge asymmetry is almost entirely model- 
independent for spin-0 and spin-1 bosons. For example, 
assuming the narrow width approximation for a Z', we 
find that Ace ^ 3/4 z{l + l/3 z'^)-l/2, for some value z 
that separates the centre and edge regions of z* , indepen- 
dent of the couplings to fermions. Spin-2 KK gravitons 
have contributions from gg and qq processes that have 
slightly different angular distributions, and the Ace de- 
pends on the weighted contribution of each. The spe- 
cific model will have an effect on the expected statisti- 
cal uncertainties, but this should not be significant to 
the measurement due to the low backgrounds associated 
with leptonic final states. Thus, with limited statistics, 
an Ace measurement could have an advantage over a fit 
to the angular distribution. 

For a hadron collider, the centre-of-mass angle of the 
outgoing fermion is not directly measurable on an event 
by event basis due to the unknown values of the par- 
ton momentum fractions. However, there exists a direct 
mapping between the cm. angular distribution and the 
difference in pseudorapidity of the final state lepton and 
anti-lepton, Ary. Furthermore, it is straightforward to 
show that At^ is a Lorentz invariant quantity, so that 
measuring this quantity in the lab frame is equivalent to 
measuring it in the centre-of-mass frame: 



A?7;afc = A?7* In 



1 



1 — cos ( 



(1) 



The normalized Ar]iab distributions for spin-0, -1, -2 res- 
onances are shown in Fig. [2l where it is clear that they 
are distinct from one another. One can therefore con- 
struct a new centre-edge asymmetry using the lab frame 
A77iab distribution in place of the cm. frame angular 
distribution. 

Using the mapping given by Eq. ([T]), we define the 
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FIG. 2; Normalized Arj distribution including detector ac- 
ceptance cuts (|77;| < 2.5, pti > 20 GeV) and only including 
events within \Mr - Mj+j- | < 0.5 AAf. These cuts reduce 
the number of measurable events with large values of |Ar;|. 
R = D (black), Z' (dark grey/orange), G (light grey/yellow), 
where only one spin-1 distribution is shown due to the model 
independent nature of the spin-1 measurement. 
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Following Osland, et al. [2J]j we take z = 0.5, which 
they find to be the "optimal" value, and translates to 
n = A77 = 1.099. Experimentally, detector acceptance 
constrains the pseudorapidity of each fcrmion to < 
2.5, which limits IA77I < 5 for the distribution. 

As is common in the literature, we take an R-parity 
violating 1/ with AA' = (0.05)^ as an example of a spin- 
resonance [l^l, and an RS graviton with c = 0.1 
for a spin-2 resonance Q. Current experimental limits 
and studies on model parameters are given in Ref. @ 
and (HI. For the Z' case we explored the Eq models ("0, 
X and T]) 1|, the Left-Right Symmetric model (LRSM, 
9R = -£l) [391], both the Littlest Higgs {LBM^an 9h = 
1.0) [4^ and Simplest Little Higgs (SLHM) [4l| models, 
and the Sequential Standard Model (SSM). 

The spin-0 model, spin-2 and some Z' models we study 
predict narrow resonances, such that including events 
within several widths of the peak will be impossible in 
practice due to detector resolution effects smearing the 
Breit-Wigner distribution. Instead, we examine events 
within one dilepton invariant mass bin as defined in the 
ATLAS TDR [Ifi], using AM = 42.9 GeV for the 1.5 TeV 
resonance as in Ref. [23] • A more precise measurement 
could be obtained by including events from a wider invari- 
ant mass window, if a broader peak were to be observed. 



TABLE I: Ace values with corresponding statistical uncer- 
tainties for 100 fb~^ integrated luminosity, pT, > 20 GeV, 
\rji\ < 2.5, within one bin AM,+;-. = 42.9 GeV and Mr = 
1.5 TeV. Also shown are the expected number of total events 
for each model assuming 100 fb^^ integrated luminosity. 



Model 


Ace ± 5 Ace 


N Events 


^6 X 
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00 / 




-0.095 ± 0.022 


2223 


Ed n 


-0.092 ± 0.021 


2480 


LR Symmetric 


-0.099 ± 0.018 


3350 


Sequential SM 


-0.097 ± 0.016 


4162 


Littlest Higgs 


-0.095 ± 0.001 


6217 


Simplest Little Higgs 


-0.094 ± 0.017 


3542 


RS Graviton 


4-0.228 ± 0.011 


8208 


R-parity violating D 


-^0.055 ± 0.066 


251 


In Table U we show the 


expected centre-edj 


;e asymme- 



try for a spin-0, spin-1 and spin-2 resonance, analogous 
to the study performed by Dvergsnes, et al. [23], assum- 
ing muon final states with 96% detection efficiency p^ . 
From Table [T] one sees that if a Z' were observed, a G 
or v could be ruled out. Likewise, an Aqe measurement 
would strongly discriminate against the Z' or v hypoth- 
esis if a G were observed. However, the Z' and G hy- 
pothesis could only be ruled out at approximately 2.5cr if 
a u signal was observed. The primary limitation in dis- 
tinguishing between the different possibilities is the low 
statistics for v production, as shown in the table, which 
is due to the tight constraints on the allowed values of 
its couplings. Other hypothetical spin-G resonances may 
not be as tightly constrained and could therefore be dis- 
tinguished from a Z' or G with higher statistical signifi- 
cance. 



IV. FORWARD-BACKWARD ASYMMETRY 

The forward-backward asymmetry is a well-established 
measurement for distinguishing between models of Z"s 
[1^ . For pp collisions at the Tevatron, the proton direc- 
tion provides an obvious choice to define the "forward" 
direction. The choice of forward direction at the LHC is 
more subtle, and is conventionally defined as the direc- 
tion of the Z' rapidity, Yz' ■ The Z' rapidity is chosen 
because the parton distribution functions for the valence 
quarks peak at a higher momentum fraction than those of 
the anti-quarks, so the system has a higher probability 
of being boosted in the quark direction. This observa- 
tion is statistical in nature and is more likely to hold 
true for larger values of ] Yz' \ ■ For smaller values of ] Yz' \ , 
the momentum fractions of the quark and anti-quark are 
generally closer in magnitude, so that using \Yzi\ in the 
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FIG. 3: Afb as a function of the Z' rapidity following Eq. (|3]) 
except that these results are not integrated over rapidity. 
From top to bottom, the models are LHM, LRSM, SLHM, 
SSM, Eeip, Eev, E^x- 

low rapidity region is less likely to correctly identify the 
quark direction. 

A simpler method of defining a "forward" or "back- 
ward" event uses pseudorapidity. As before, we define 
the quark direction to be that of the higher momentum 
parton, or equivalently the direction of the Z' rapidity. 
One can then show that a "forward" event is one in which 
1*7/ 1 > IVfl i'^ the lab frame, and vice-versa for a "back- 
ward" event. Using these definitions for forward and 
backward, one can define the forward-backward asym- 
metry: 



Afb = 



[F{y) - B{y)] dy 
[Fiy) + Biy)] dy 



(3) 



where F(y) is the number of forward events and B{y) is 
the number of backward events for a given y, the Z' ra- 
pidity (i.e.: Yz'). The F{y)—B{y) distribution under this 
definition is clearly shown in Fig. [3] to be symmetric in 
Z' rapidity. This method of finding Apg has the advan- 
tage of being very straightforward and clean. It simply 
relies on counting events with |?7y| > \rijr\ and those with 
\Vf\ < I '7/ 1- We note that a related technique is employed 
by the CDF collaboration [i^l for the Z° ApB inpp col- 
lisions at the Tevatron. However, the natural choice of 
the quark direction in pp collisions at the Tevatron in 
contrast to pp collisions at the LHC results in important 
differences between the methods. 

As in the conventional method for finding ApB, for 
small values of there is a higher probability to 

wrongly assume that the quark is the parton with the 
higher momentum fraction. This results in incorrectly 
assigning the forward or backward direction and gives 
a small "wrong" contribution to the Afb measurement. 
For this reason, it has been suggested that the central re- 
gion, \Yz' \ < Ymin, be excluded in the measurement [1^. 
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FIG. 4: Afb off-peak versus on-peak for a variety of models, 
including detector acceptance limits and kinematic cuts as 
previously listed. Standard Model measurement determined 
from the standard model Drell-Yan cross section, with on- 
peak events within |M;+;- — Af^/ 1 < 100 GeV and off-peak 
events within 2/3Mz/ < M,+;- < Mz' - 300 GeV to include 
large enough statistics. 



However, the coupling dependency can still be deter- 
mined without this constraint on \Yzi\ [l9l |. 

Another consideration for excluding the central region 
is that the number of events that remain after subtract- 
ing F — B \s small, as shown in Fig. (H while the total 
number of events in this region is large. Excluding the 
events in the central region would increase the magnitude 
oi ApB^ potentially making models more distinguishable. 
However, we found that increasing Ymin resulted in an in- 
crease in the relative uncertainty. We therefore conclude 
that little is gained by excluding events with small Y^/, 
and suggest that the whole rapidity region be included to 
decrease uncertainty and further simplify the ApB mea- 
surement. 

Using this method, we calculate Afb for the i^g mod- 
els {tpjXjV) the Left Rig ht Symmetric model [391, 
the Littlest Higgs model ^0|, the Simplest Little Higgs 
model [iH , and the Sequential Standard model. The on- 
peak versus off-peak Afb are shown in Fig. ID, where on- 
peak includes events which satisfy |M;+;- — Mz' \ < 3Tz' 
and off-peak includes events which satisfy 2/3Mz> < 
Mi+i- < Mz' — Sr ^/ , similar to the cuts used by Petriello 
and Quackenbush |20 |. 

We conclude with the important observation that it 
might also be possible to include some events in the for- 
ward regions of the calorimeter (FCAL) using this tech- 
nique [43[. While a muon signature appears as missing 
Et in the FCAL, it may be possible to distinguish an 
electron from a jet in the FCAL due to the differences in 
the showering. The signal would require triggering off of 
a single, high px electron in the jryj < 2.5 region, with an 
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electron-jet in the FCAL. Determining the charge sign 
of the single electron would distinguish whether this is 
a "forward" or "backward" tagged event. It is not clear 
what the signal efficiency of this method is, as reducible 
backgrounds include W + j and others that might have 
low rejection rates. Extending the rapidity range has 
the potential of increasing the statistics and remains an 
interesting possibility for further study. 

V. SUMMARY 

In this paper we described an approach for discriminat- 
ing between various spin hypotheses for a newly discov- 
ered s-channel resonance at the LHC using a centre-edge 
symmetry, Ace, that is based on the difference of the 
rapidities of final state fermions. We also described a 
simple way to measure the forward-backward asymme- 



try, ApB, using the properties of pseudorapidity. Both 
of these measurements have an advantage over previous 
approaches as they rely solely on the measurement of 
pseudorapidity, a fairly basic quantity. The new mea- 
surements require simple counting and should propagate 
fewer errors than previous approaches that rely on boost- 
ing the four momentum into the centre-of-mass frame in 
order to perform the analysis. Our approaches repro- 
duce the results found in other analyses but via a more 
straightforward analysis. 
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